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Abstract 

The standing MHD modes in a zero-/3 cyhndrical magnetic flux tube modelled as a 
straight core surrounded by a magnetically twisted annulus, both embedded in a straight 
ambient external field is considered. The dispersion relation for the fast MHD waves is 
derived and solved numerically to obtain the frequencies of both the kink (m = 1), and 
fluting (m — 2, 3) waves. Damping rates due to both viscous and resistive dissipations in 
presence of the twisted magnetic field is derived and solved numerically for both the kink 
and fluting waves. 
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1 Introduction 

Solar corona is a highly structure of magnetic flux tubes namely coronal loops. Transverse 
oscillations of coronal loops were first identified by Aschwanden et al. (1999) and Nakariakov 
et al. (1999) using the observations of TRACE. Edwin & Roberts (1983) elaborated on the 
dispersion relation for a magnetic cylinder embedded in a magnetic environment typical of that 
of the solar photosphere and corona. They found that the existence of inhomogeneities in the 
form of structuring of the magnetic field enables loops to act as wave guides for a variety of 
different modes. Karami, Nasiri & Sobouti (2002) used the model of Edwin & Roberts (1983) 
but without limiting it to slender tubes. They solved numerically the dispersion relation for 
each mode in its full generality. They obtained that in presence of weak viscous and ohmic 
dissipations, the damping rate is inversely proportional to the Reynolds and Lundquist numbers, 
R and S, respectively. 

An additional features of the flux tube is that of twist. Bennett, Roberts & Narain (1999) 
examined the influence of magnetic twist on the modes of oscillations of a magnetic flux tube. 
They found that twist introduces an infinite band of body modes. Klimchuk, Antiochos & 
Norton (2000) introduced twist to resolve the internal structure on an individual loop embedded 
within a much larger dipole configuration. Mikhalyaev & Solov'ev (2005) investigated the MHD 
waves in a double magnetic flux tube embedded in a uniform external magnetic field. The tube 
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consists of a dense hot cylindrical cord surrounded by a co-axial shell. They found two slow and 
two fast magnetosonic modes can exist in the thin double tube. 

Verwichte et al. (2004), using the observations of TRACE, detected the multimode oscil- 
lations for the first time. They found that two loops are oscillating in both the fundamental 
and the first-overtone standing fast kink modes. According to the theory of MHD waves, for 
uniform loops the period ratio P1/2P2 of the fundamental mode and its first overtone is exactly 
1. But the ratios found by Verwichte et al. (2004) are 0.91±0.04 and 0.79±0.03 and thus clearly 
differ from 1. This may be caused by different factors such as the effects of curvature, see e.g. 
Van Doorsselaere et al. (2004), leakage, see De Pontieu et al. (2001), density stratification in 
the loops, see e.g. Andries et al. (2005), Erdelyi & Verth (2007), Karami k Asvar (2007) and 
magnetic twist, see Erdelyi & Fedun (2006) and Erdelyi &; Carter (2006). 

Erdelyi & Fedun (2006), studied the wave propagation in twisted cylindrical magnetic flux 
tube embedded in an incompressible but also magnetically twisted plasma. They found that 
magnetic twist will increase, in general, the periods of waves approximately by a few percent 
when compare to their untwisted counterparts. Erdelyi &: Carter (2006) used the model of 
Mikhalyacv &: Solov'ev (2005) but for a fully magnetically twisted configuration consisting of 
a core, annulus and external region. They investigated their analysis by considering magnetic 
twist just in the annulus, the internal and external regions having straight magnetic field. Two 
modes of oscillations occurred in this configurations; surface and hybrid modes. They found 
that when the magnetic twist is increase the hybrid modes cover a wide range of phase speeds, 
centered around the annulus, longitudinal Alfven speed for the sausage modes. 

Carter k. Erdelyi (2007) investigated the oscillations of a magnetic flux tube configuration 
consisting of a core, annulus and external region each with straight distinct magnetic field in 
an incompressible medium. They found that there are two surface modes arising for both the 
sausage and kink modes for the annulus-core model where the monolithic tube has solely one 
surface mode for the incompressible case. Also they showed that the existence and width of an 
annulus layer has an effect on the phase speeds and periods. Carter &: Erdelyi (2008) used the 
model introduced by Erdelyi & Carter (2006) to include the kink modes. They found for the set 
of kink body modes, the twist increase the phase speeds of the modes. Also they showed that 
there are two surface modes for the twisted shell configuration, one due to each surface, where 
one mode is trapped by the inner tube, the other by the annulus itself. 

In the present work, our aim is to investigate the effects of the twisted magnetic field on 
oscillations and damping of standing MHD waves in the cold coronal loops observed by Verwichte 
et al. (2004) deduced from the TRACE data. This paper is organized as follows. In Section 
2 we use the model introduced by Erdelyi & Carter (2006) to derive the equations of motion, 
introduce the relevant boundary conditions and obtain the dispersion relation. In Section 3 
we discuss resistive and viscous dissipations to calculate contributions of the different modes to 
heating of the coronal loops. In Section 4 we give numerical results. Section 5 is devoted to 
conclusions. 



2 Equations of Motion 

The linearized MHD equations for a zero-/? incompressible plasma are 

^ = 7^{(V X (5B) X B + (V X B) X (5B} + ^V^Jv, (1) 
at Airp p 

^ = Vx(«vxB) + _V^«B, (2) 
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(3) 



V-6v = 0, (4) 

where 5v, (5B and Pt are the Eulerian perturbations in the velocity, magnetic field and magnetic 

pressure, respectively; p, a, i] and c are the mass density, the electrical conductivity, the viscosity 
and the speed of light, respectively. Note that Eq. (4) satisfies the incompressibility condition. 
The simplifying assumptions are: 

• p is constant over the loop; 

• for a zero-/? loop, gas pressure is negligible; 

• from Erdelyi &; Carter (2006), background magnetic field is assumed to be 



B = I 



Bi = {0,0, Bi), r<a, 

Bo = {0,Aor,Bo), a < r < R, 

Be = {0,0, Be), r>R, 



where Bi, Bq, Be are constant and a, R are radii of the core and tube, respectively; 

• tube geometry is a circular with cylindrical coordinates, {r,4>,z); 

• there is no initial steady flow over the tube; 

• viscous and resistive coefficients, r] and a respectively, are constants; 

• t-, 0- and z- dependence for any of the components (5v and 5B is exp {i{m(j) + kzZ — uot)}. 
Where kz = Itt/L, L is length of the tube, and Z = (1,2, • • •), m = (0, 1,2, • • •) are the 
longitudinal and azimuthal mode numbers, respectively. 

We will further assume that the dissipative terms in Eqs. (1) and (2) are much smaller. We 
will first solve the problem without these terms and re-introduce them later as small correc- 
tions in calculating contributions of the different modes to heating of the corona. Taking time 
derivative of Eq. (2) and substituting for dS\/dt from Eq. (1), the resulting equation yields to 
Bessel's equation 

dr^ r dr 

where 

2 



-{^ + ml)PT = 0, (5) 



^2 



and 



^Ao = + kzBo). (7) 

Equation (5) is same as the result exactly derived by Erdelyi & Carter (2006). Note that 
subscripts (which correspond to annulus) are replaced by i, e corresponding to the internal 
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and external regions, respectively. Since for the internal and external regions Ai=Ag=0, hence 
m,? = = A;2 > and ua, = " 



Solutions of Eq. (5) are: 

Pt = alm{k,r), (8) 

for the interior region (r < a), 

^ r pirnimor) +-fKrnimor), > 0, 

I pJminor) +'yYm{nor), = -m^ > 0, ^ 

for the annulus region [a < r < R) and 

Pt = 5Xm(A;.r), (10) 

for the exterior region {r > R). Where {Jm-,Ym) and (J^, K^^) are the Bessel and modified Bessel 
functions of the first and second kind, respectively. The coefficients a, /3, 7 and 5 are determined 
by the boundary conditions. From both Karami, Nasiri & Sobouti (2002) and Erdelyi & Carter 
(2006), the necessary boundary conditions are that: at the boundaries r = a and r = R, both 
the Lagrangian magnetic pressure and dvr should be continuous. These conditions yield to the 
dispersion relations for surface, mg > 0, and hybrid, mg < 0, modes which are same as the 
results obtained by Erdelyi & Carter (2006) in Eqs. (28a) to (28b), respectively. Note that 
numerical solution of the dispersion relation yields to eigenfrequencies, which are characterized 
by a trio of wave numbers (n, m, I) that actually count the number of nodes or antinodes along 
r, (f), and z directions, respectively. 



3 Dissipative Processes 

Since the discovery of the hot solar corona about 66 years ago, different theories of coronal 
heating have been put forward and debated. For instance, Nakariakov et al. (1999) reported 
the detection of spatial oscillations in five coronal loops with periods ranging from 258 to 320 s. 
The decay time was (14.5 it 2.7) minutes for an oscillation of (3.9 it 0.13) millihcrtz. Also Wang 
k Solanki (2004) described a loop oscillation observed on 17 April 2002 by TRACE in 1951. 
They interpreted the observed loop motion as a vertical oscillation, with a period of 3.9 minutes 
and a decay time of 11.9 minutes. All these observations indicate strong dissipation of the wave 
energy that may be the cause of coronal heating. 

Following Karami & Asvar (2007) , the finite conductivity and viscosity of plasma causes an 
exponential time decay of disturbances. Hence for weak dissipations one may assume 

j-gdissipative _ (5B(r)e~*^*'^"'"")* , 

^^dissipative ^ 5v(r)e-(^'^+")*, (11) 

where lo, d'B, and 5v on the right hand side are the solutions of Eqs. (1), (2) in the absence of 
dissipations. Substituting Eq. (11) in Eqs. (1) and (2), canceling out the non dissipative terms, 
and keeping only the first order terms in a, (? /A-ku and r] gives 

2iawJB = - — I [(B • V)(^v] - [((5v • V)B] [ 

+ ^{(B-V)V2(5v-(V^(5v-V)b}, (12) 
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where 



Sv = - 



UJ 



UJz 



-{v(B-5B)-((5B-V)b}. 



(13) 



Rewriting Eq. (12) for either the transverse or the z-component and substituting for ah 
quantities in terms of SB^ gives 



a 



c 



Ana po 



Po V27r(a;'' 



where va. 



Bi 



the Lundquist number S 



( 47rg-R^ 



1 1 



(14) 



27rR 



is the ratio of the resistive time 



27rR 



is the ratio of 



scale to the Alfven crossing time and the Reynolds number TZ 

the viscous time scale to the Alfven crossing time. Equation (14) shows that when the twist is 
absent, Aq = 0, the damping rate is vanished. Whereas for compressible plasma it is not zero. 
See Karami, Nasiri & Sobouti (2002). 



4 Numerical Results 

As typical parameters for a coronal loop, we assume L = 109 X 10^ km, R/L = 0.01, pe/pi = 0.1, 
po/pj = 0.5, pi = 2x 10-^^ gr cm-3, Be/Bi = 1, Bo/Bi = 1, 5j = 100 G. For such a loop one 
finds VA, = 2000 km s'^, loa^ := ^ = 1-835 rad s"^. We use S = 10^ and 7^ = 560, given by 
Ofman et al. (1994). 

The effects of twisted magnetic field on both the frequencies oj and damping rates a are 
calculated by numerical solution of the dispersion relation, i.e. Eqs. (28a)- (28b) given by Erdelyi 

& Carter (2006), and Eq. (14), respectively. The results are displayed in Figs. 1 to 7. Figures 
1 to 6 show the frequencies and damping rates of the fundamental and first-overtone / = 1,2 
kink (m = 1) and fluting (m = 2,3) surface modes versus the twist parameter, B^/Bz := 
and for different relative core width a/R = (0.65,0.8,0.9). Figures 1 to 6 reveal that: i) For 
a given a/R, both frequencies and their corresponding damping rates increase when the twist 
parameter increases. The result of lo is in good agreement with the that obtained by Carter 
k, Erdelyi (2008). Note that there are two surface modes labelled by (n = 1,2) which are in 
accordance with Carter k, Erdelyi (2008). Here we only show the first one (n = 1) in the figures, 
because the second one {n = 2) does not show itself in all selected twists, ii) For a given m and 
a/R, when the longitudinal mode number, increases, both the frequencies and damping rates 
increase, iii) For a given /, a/R and B^p/Bz, when the azimuthal mode number, m, increases, 
the frequencies and damping rates increase and decrease, respectively. 

Figure 7 presents the frequencies and damping rates of the kink {m = 1) surface modes with 
I = 100 versus the twist parameter. Figure 7 shows that for I = 100, the damping becomes 
stronger and the ratio u/a decreases two order of magnitude compared with I = 1, 2. See again 
Figs. 1 to 2. 

Here in our calculations, the sausage modes (m = 0) are absent. Because following Edwin 
& Roberts (1983) and Karami, Nasiri & Sobouti (2002), the sausage modes have a lower lon- 
gitudinal cutoff and they are only expected in fat and dense loops. For instance, according to 
Aschwanden (2005) for typical active region loops which have a density contrast in the order of 
Pe/Pi ~ 0.1 — 0.5, would be required to have width-to-length ratios of L/(2R) 1 — 2. 
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The period ratio P1/2P2 of the fundamental and first-overtone, I = 1,2 modes of both the 
kink (m = 1), and fluting (m = 2,3) surface waves versus the twist parameter plotted in Figs. 
8 to 10, respectively. Figures 8 to 10 show that: i) For a given relative core width, the period 
ratio P1/2P2 decreases when the twist parameter increases. For a/i?=0.65, for instance, P1/2P2 
decrease from 1 (for untwisted loop) and approaches below 0.95, 0.88 and 0.82 for m = (1,2,3), 
respectively, with increasing the twist parameter, ii) For a given twist parameter, the period 
ratio P1/2P2 increases when the relative core width increases. Figure 8 clears that for kink 
modes (m = 1) with 5^/5^=0.0065 and a/R = 0.65, the ratio P1/2P2 is 0.941. This is in good 
agreement with the period ratio observed by Verwichte et al. (2004), 0.91±0.04 deduced from 
the observations of TRACE. See also McEwan, Diaz & Roberts (2008). 

Figure 11 displays the frequency band width. Aw, including infinite set of the fundamental 
kink (m = 1) hybrid modes versus the twist parameter and for different relative core width. 
Figure 11 presents that: i) For a given twist parameter, Ao; increases when the relative core 
width decreases, ii) For a given relative core width, Ao; increases when the twist parameter 
increases. This is in good agreement with the result obtained by Carter &; Erdelyi (2008). 

5 Conclusions 

Oscillations and damping of standing fast MHD surface and hybrid waves in coronal loops in 
presence of twisted magnetic field is studied. To do this, a typical coronal loop is considered 
as a straight pressureless cylindrical incompressible flux tube with magnetic twist just in the 
annulus and straight magnetic field in the internal and external regions. The linearized MHD 
equations, when the dissipation is absent, are reduced to a Bessel's equation for the perturbed 
magnetic pressure. The dispersion relation is obtained and solved numerically for obtaining 
the frequencies of both the kink and fluting modes. The damping rates of oscillations due to 
the resistive and viscous dissipation in presence of the magnetic twist is obtained and solved 
numerically. Our numerical results show that: 

i) For a given relative core width, frequencies and damping rates of both the kink (m = 1) 
and fluting (m = 2, 3) surface waves increase when the twist parameter increases. 

ii) The period ratio P\j2P2, for both the kink (m = 1) and fluting (m = 2, 3) surface modes 
are lower than 1 (for untwisted loop) in presence of the twisted magnetic field. The result of 
P1/2P2 for kink modes is in accordance with the TRACE observations. 

iii) Frequency band width of the fundamental kink (m = 1) hybrid modes increase when the 
twist parameter increases. 
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Kink Surface Modes (m=1, 1=1) 




xio ■ 




Figure 1: Frequency of the fundamental kink (m = 1) surface mode and its damping rate versus 
the twist parameter, B^j^/Bz, for different relative core width a/R = 0.65 (solid), 0.8 (dashed) 
and 0.9 (dash-dotted). The loop parameters are: L = 109 x 10^ km, R/L = 0.01, Pe/ Pi = 0.1, 
pa/pi = 0.5, Pi = 2 X lO"^'' gr cm^^ Be/Bi = 1, Bo/Bi = 1, Bi = 100 G, S = 10^ and 
TZ = 560. Both frequencies and damping rates are in units of the interior Alfven frequency, 
u>Ai = 1.835 rad s"^. 
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Kink Surface Modes (m=1, 1=2) 




Figure 2: Same as Fig. 1, for the first-overtone kink (m = 1) surface modes. 



Fluting Surface Modes (m=2, 1=1) 




Figure 3: Same as Fig. 1, for the fundamental fluting (m = 2) surface modes. 
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Fluting Surface Modes (m=2, 1=2) 
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Figure 4: Same as Fig. 1, for the first-overtone fiuting (m = 2) surface modes. 
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Fluting Surface Modes (m=3, 1=1) 
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Figure 5: Same as Fig. 1, for the fundamental fluting (m = 3) surface modes. 
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Fluting Surface Modes (m=3, 1=2) 
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Figure 6: Same as Fig. 1, for the first-overtone fiuting (m = 3) surface modes. 
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Figure 7: Same as Fig. 1, for the kink (m = 1) surface modes with I = 100. 



11 



Kink Surface Modes (m=1) 
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Figure 8: The period ratio P1/2P2 of the fundamental and its first-overtone kink (m = 1) 

surface modes versus the twist parameter for different relative core width a/R = 0.65 (solid), 
0.8 (dashed) and 0.9 (dash-dotted). Auxiliary parameters as in Fig. 1. 



Fluting Surface Modes (m=2) 




♦ z 

Figure 9: Same as Fig. 8, for fluting (m = 2) surface modes. 
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Fluting Surface Modes (m=3) 
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Figure 10: Same as Fig. 8, for fluting (m = 3) surface modes. 



Kink Hybrid l\^odes (m=1 , 1=1 ) 




Figure 11: Frequency band width of tlie fundamental kink (m = 1) hybrid modes versus the 
twist parameter for different relative core width a/R = 0.65 (solid), 0.8 (dashed) and 0.9 (dash- 
dotted). Auxiliary parameters as in Fig. 1. 
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